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Abstract 

We study non-perturbatively and from first principles the thermodynamics of vortices 
in 3d U(l) gauge+Higgs theory, or the Ginzburg-Landau model, which has frequently 
been used as a model for cosmological topological defect formation. We discretize 
the system and introduce a gauge-invariant definition of a vortex passing through a 
loop on the lattice. We then study with Monte Carlo simulations the total vortex 
density, extract the physically meaningful part thereof, and demonstrate that it has 
a well-defined continuum limit. The total vortex density behaves as a pseudo order 
parameter, having a discontinuity in the regime of first order transitions and behaving 
continuously in the regime of second order transitions. Finally, we discuss further 
gauge-invariant observables to be measured. 
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1 Introduction 



Vortices play a significant role from the low temperatures of liquid crystals |]J , super- 
fluids @ and high-T c superconductors || , to the relativistic temperatures of the Early 
Universe [|J. In low temperature systems, vortices can be directly observed |T], [J; in 
cosmology, one has studied their effects on the inhomogeneities leading to structure 
formation ||. Consequently, vortices have been a subject of immense interest during 
the last few years. Nevertheless, some important questions, related in particular to 
non-perturbative studies of vortices in gauge theories, remain poorly understood. 

Among the most fundamental principles of Nature appear to be gauge invariance and 
spontaneous symmetry breaking, and already the simplest theory with these properties, 
the locally U(l) symmetric gauge + scalar quantum field theory or the Ginzburg-Landau 
(GL) model, contains vortices. The GL model does describe real physics in liquid 
crystals and superconductors 0, while in cosmology it is to be viewed as a simple toy 
model. The phase structure of the GL model is non-trivial: in the type I regime, there 
is a first order transition f8[, whereas in the type II regime, the transition is assumed 



to be of the second order Hl-112 



One of the mentioned open questions arises immediately when one realizes that the 
type II regime of the GL model is completely non-perturbative: perturbation theory 
does not describe the transition at all |13| |. The only known systematic and controllable 
method for studying this regime are lattice Monte Carlo simulations ||13||-||18||. Yet as to 
date, to our knowledge, the vortex density has not been studied in detail on the lattice 
even in thermodynamical equilibrium. The purpose of this paper is (a) to provide a 
gauge-invariant formulation for studying vortices on the lattice, (b) to measure the 
vortex density both in type I and type II regimes, and (c) to extrapolate the results 
to the infinite volume and continuum limits. The length distribution of vortices will 



be studied in a future publication [fL9|| . In a U(l) scalar field theory without gauge 
symmetry, the thermodynamics of vortices has previously been addressed in |[20|| . 

Let us stress that considering the thermodynamics of vortices is certainly only a 
starting point. Ultimately one is interested in the real-time scaling properties of vortex 



networks created in a non-equilibrium situation (see, e.g., |21[|-[p3|l). However, the 
thermodynamical equilibrium situation provides the initial conditions for such non- 
equilibrium processes. Furthermore, it is clear that non-equilibrium physics cannot be 
understood in quantitative detail before the equilibrium limit is under control. 



2 The theory in the continuum and on the lattice 

Let us start by defining the theory. The continuum theory is defined by the functional 
integral 

Z = |pA i P0exp[-S'(A,0)], (1) 
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Figure 1: The phase diagram of the GL model, together with the main simulation 
points (marked with the diamonds). 
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where = diAj — djA t and Di = di + ie%Ai. The theory is invariant under the gauge 
transformations 



3 i0(x) 



(x), A,(x) -> A 4 (x) - c^(x)/e 3 . 



(3) 



Writing 0(x) = v{x) exp[z7(x)], the first of these can be rewritten as 7(x) — > [7(x) + 
^(x)],,-, where [X] n = X+27rn such that [X] n G (— 7r, tt] and we have chosen to represent 
the phase of by a number in this interval. The theory in eq. (H) is parameterized by 
the scale and by the two dimensionless ratios 
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where Tn\{^i) is the mass parameter in the MS dimensional regularization scheme in 
3 — 2e dimensions. Expressions for x and y in terms of the original physical parameters 
of both 4d high temperature scalar+fermion electrodynamics and 3d low-T c supercon- 



ductivity, have been discussed in [pl|] where we refer to for more details. Here we just 
study the theory as a function of x, y. The phase diagram is shown in Fig. [TJ. 

As is well known, the classical counterpart of the theory in eq. (0) admits vortex, or 
string solutions, in the broken (superconducting) phase. In other words, the classical 
equations of motion have solutions in which the symmetry is restored at the core, but 
is broken far away from the core p4| . The existence of a vortex inside a loop C can be 



identified by computing the line integral 

l£dx.V 7 (x), 



n c 



(5) 
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where the winding number nc is an integer, and a non-zero nc signals a vortex inside 
the loop. It can be seen from eq. ([|) that nc is gauge-invariant for single- valued gauge 
transformations 0(x). Thus vortices are physical objects, which can be generated 
in a phase transition or by applying non-trivial boundary conditions, and they can 
be observed in the superconducting phase. However, in the equilibrium state they 
are also created and destroyed by thermal fluctuations, and (|^c|) gives their density 
both in the broken and in the symmetric phase. Moreover, it is expected that the 
behavior of vortices is qualitatively different in the two phases. As the dynamics is 
non-perturbative, it is then clear that it is very important to be able to study vortices 
with lattice simulations. 

It should be noted here that one often assumes that in the regime of large x, the 
gauge fields are not essential and one can approximate the theory in eq. (§) by the 3d 
XY-model with a global U(l) symmetry, or its dual version (see, e.g., |L3], whose 
fundamental objects are the vortices. While these approximations simplify the problem 
significantly, their validity is uncertain for finite x. Thus we consider it essential to 
approach the problem directly with the original theory in eq. (12). 

To allow for lattice simulations, the theory in eq. (Q) has to be discretized. As 
usual, we introduce the link field Z7j(x) = exp[mesA(x)] = exp[ic^(x)]. Rescaling the 
continuum scalar field to a dimensionless lattice field by <p*(f) —>■ /?#0*</>/2a, the lattice 
action becomes 

s = PgT, ^'( x )-^E Re ^( x )^( x )^( x + 

X,i<j Z X,i 

+ £ 0* (x)0(x) + (3 R £ [0*(x)0(x) - l] 2 , (6) 

X X 

where ^(x) = CKj(x) + atj(x + 1) — otj(x + j) — ct/(x). We use here the non-compact 
formulation for the gauge fields. The gauge transformation properties in eq. (^) go over 
into 

7 (x) - [ 7 (x) + 0(x)] 7r , 

ai(x) -> a<(x) + 0(x)-0(x + i). (7) 

Discretization can be viewed as a different regularization scheme, and in order to de- 
scribe the same continuum physics as in eq. (0), one has to make a 2-loop computation 
relating the counterterms [p8[| . As a result, the lattice couplings Pg,Ph,Pr are deter- 
mined from 





, ; , •■ u , = _ 3.1759115(1 + 2x)f3 G 
-4 + 8x - 8a; 2 )(log6/3 G + 0.09) -1.1 + 4.6x 

167T 2 
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Thus for a given continuum theory depending on one scale e| and the two dimensionless 
parameters y, x, the use of a lattice introduces a regulator scale a, and eqs. (|6])-(|9]) 
specify, up to terms of order e^a, the corresponding lattice action. Note that the 
simulations in [^5| correspond to (3r — > oo, whereas according to eq. /3r — > in 
the continuum limit for any finite x {(3h — > 1/3). 



3 Gauge- invariant vortices 

Consider now vortices. The naive discretization of eq. @ gives the standard algorithm 
used in scalar theories without gauge fields. For each loop C one would define the 
winding number he of the phase 7 of the scalar field. However, any 7(x) can be changed 
arbitrarily with a gauge transformation, see eq. ([T|). Thus the 7(x)'s are essentially 
random numbers, and he does not contain any real information about the dynamics of 
the system. Indeed, the result for a loop around a single plaquette, iVi x i = (|ni x i(x)|), 
would equal 1/3 in the case of completely uncorrelated field values, and this is what we 
measure from lattice simulations for AT lxl , irrespective of the parameters of the theory 
(to be more precise, we always get iVi x i = 0.32(1) . . .0.33(1)). Thus the quantity he 
has to be rejected. 

One solution sometimes used in the literature would be to fix the gauge, which makes 
h c non-trivial. However, its value depends crucially on the gauge chosen, and it is even 
possible to choose a gauge in which he always vanishes. Therefore, we believe that 
it is important to use an explicitly gauge-invariant definition, in order to be able to 
interpret the results correctly. 

Fortunately, the problem with he is not one of principle, and a satisfactory definition 
can be given. For each positively directed link I — (x, x + i) let us define 

Y {x , x +t) = + 7(x + i) - 7 (x)] ff - Oi(x). (10) 

For links with negative direction the sign of Yi is changed: V( XjX _i) = — ^( x -?,x)- Then, 
for each closed loop C, we can define 

Y c = Y. Y i = ^ n c- (11) 

zee 

This definition has four main properties: 

(a) For any field configuration and any loop C, nc € 2. 

(b) Directly from eqs. (0), one can see that the part of Yj in the square brackets is 
gauge-invariant. The term — a t -(x) is not, but when summed over a closed loop into 
Ye, the gauge dependence cancels. Hence Ye is gauge-invariant. 

(c) Since Ye is gauge-invariant, one can always tune the gauge used in the evaluation 
of the Yi's such that the fields appearing are perturbatively small. But then, in the 



4 



continuum limit, at = e^aAi goes to zero and one gets the correct continuum limit 
containing only the phase of 0. 

(d) The quantity Yc is additive: if there is a loop C consisting of the loops A, B, then 
Yc = Ya + Yb- This is what one would require for counting the number of vortices 
going through loops of different sizes. This also implies that vortex lines cannot end, 
and therefore they form closed vortex loops. Note that we use a non-compact gauge 
field so that there are no monopoles. 

Based on these properties, eq. (|TTD provides a valid formulation for counting vortices 
in the locally symmetric U(l) theory |30| . This definition of the winding number 



coincides with that given in Ref. [[29j for the case /3r = oo. 



4 Simulations and results 

In order to see how the gauge-invariant definition performs in practice, we have made 
lattice Monte Carlo simulations in the GL model. We choose C = n x n = a loop 
around a plaquette of size nxnin eq. (|ITD, and measure 

N nxn = (|n nxn |)- (12) 

The quantity N nxn measures the average net number of vortices through a loop of size 
nxn, irrespective of the net direction. Keeping track of the direction would give zero: 
for symmetry reasons, (n nxn ) = 0. In practice, we average (\n nxn \) over all lattice sites 
and directions, to improve on the statistics. 

Simulations are made at two values of x: x = 0.0463 corresponds to a strongly type I 
superconductor, x = 2 to a strongly type II superconductor. For each x, values of y are 
chosen on both sides of the transition (see Fig. [I]). For each such continuum parameter 
point, several lattice spacings are chosen: [3g = l/(e|a) = 1, 2, 3, 4, 6, 8, 12. For j3c = 4 
(x = 0.0463) and fla = 1,2 (x — 2), several volumes are chosen, in order to test that 
the finite volume effects are small. The volumes thus arrived at (24 2 x48 for (3g = 4 
both at x — 0.0463, x = 2) are then scaled with (3q such that the physical volume (in 
units of l/eg) remains constant. The orders of magnitude for j3c and the volume come 
from the requirement that the physical correlation lengths, of order (0.5...2)/e| at 
x = 0.0463 and (1.5 . . . 3)/e| at x = 2 |I7 |, are much longer than the lattice spacing 



but much shorter than the extent of the whole lattice (with the exception of the photon 
in the symmetric phase). 

The results for Ni function of y are shown in Fig. The two values of x have 

a different critical point y c , see Fig. [I]. At small x the vortex density is very small in the 
broken phase but jumps discontinuously to a large value at the transition. However, 
the symmetric-phase value ~ 0.2 is much smaller than the trivial value 1/3. At large x 
the behavior is completely different. The total vortex density is large also rather deep 
in the broken phase and there is no discontinuity at the transition. 
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Figure 2: The values of iV lxl at x — 0.0463 and x = 2, as a function of y. The 
transition points are indicated by the arrows. There is a discontinuity at the first order 
transition at x = 0.0463, whereas in the regime of large x, iV lxl behaves continuously. 
Some of the points at x = 0.0463 correspond to metastable phases. 



Let us then discuss the approach to the continuum limit: a — > and hence (3g — > oo. 
The continuum extrapolation of iV lxl is shown in Fig. |^. It is seen that even though 
there is a lot of structure at a finite /?g, all the structure disappears when [3g — > oo and 
one gets a result which is independent of the parameters of the continuum theory [31]. 
In the continuum limit, the loop of size lxl (as well as a loop of any finite size nxnin 
lattice units) shrinks to a point, and it is clear that the result is merely an artifact of the 
lattice regularization, and sensitive only to UV-effects. Note that the continuum value 
Nixi ~ 0.2 is close to the "universal" value seen, e.g., in pOf (there it was obtained 
at a fixed lattice spacing but at the point of a second order phase transition, which 
corresponds precisely to the continuum limit). 

An important point to be noticed from Fig. |3| is that in the broken phase (x = 
0.0463,?/ = 0.08), the asymptotic /3g regime is obtained quite late, /3g^,8. This is 
somewhat surprising since the smallest physical correlation length at this point is 
1/mw ~ 0.4/e| [|T7|, corresponding to ~ 2.4a already at Pg — 6. Thus one would 
expect to be approaching the continuum limit earlier. The Higgs correlation length is 
much larger, l/m H ~ 2/e^ ~ 12a at f3 G — 6. 

For larger loops nxn, the qualitative behaviour is similar to that for iV lxl , although 
the numerical values are different. For large /3g, including only terms linear in a, we 
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Figure 3: The continuum extrapolation of iV lxl for different values of x,y. Even 
though there is a lot of structure at finite /?g, all results approach a single point in 
the continuum limit (3g — > oo, including those in the broken phase. Both polynomial 
(co + ci/Pg + c 2/Pg) an d logarithmic (co + ci\d.^g/Pg + C2/A3) fits are allowed at 
1//3g - ► 0, the latter ones favouring a slightly larger Co ~ 0.20. 

expect the values of N nxn to behave as 

N nxn « /(n) + Kx, ?/) + e(x, y)n]//3 G + 0(1/^). (13) 

In principle, there could be a ln/3 G -term in <i(x, y), e(x, y). For fixed n, the continuum 
value is f(n), but it does not reflect the infrared dynamics of the theory. Fitting the 
data, the functional form of f(n) is found to be consistent with a + b/n + c Inn for 
large n. We cannot conclusively determine whether the coefficient c is non-vanishing 
or not. Assuming c = 0, we get a ~ 0.33, b ~ —0.13, but if c is allowed to be non-zero, 
the absolute values of a and b are somewhat smaller. The numerical determination of 
c is difficult, since it requires large values of n, for which very large lattices are needed 
to remove the finite size effects. In any case, the real physics lies in the coefficient of 
1/Pg ( see below), for which a fit of the form in eq. ( |T3"D works very well (the confidence 
level is CL=10-90%, depending on the parameter values x, y). 

Let us then discuss loops which are of a fixed size in physical units: (c/e|) x (c/ef), 
where c is a constant. In lattice units this loop is c/3g x cPg, i-e., the size is varied as 
a function of (3q- The continuum extrapolation for the loop Np aX p G (c = 1) is shown 
in Fig. [| at the point x = 0.0463,?/ = 0.14. According to eq. (|13D , we expect that in 
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Figure 4: The continuum extrapolation of Np GX p G at x — 0.0463,?/ = 0.14. Linear, 
quadratic and logarithmic fits are shown; the logarithmic one has the best confidence 
level, CL=11%. The value of Np GX p G in the continuum limit is dominated by an 
unphysical regularization sensitive constant part, see the text. 

the continuum limit, 

N {c p G)xic( 3 G) ~ lim f(c(3 G ) + ce(x, y) + 0(c 2 ). (14) 

The first term is unphysical and corresponds to the regularization effects in eq. (|i~5|) 
below. (If there is a term ln/^ in e(x,y), then the regularization sensitive part can 
also depend on x, y, but this dependence is analytic and does not affect any phase 
transitions.) Thus the absolute value of Nf c p G ) x r c p G ^ is not physical, only its changes 
are (see Fig. ||). The physical effects, i.e. ce(x, y) + C(c 2 ), come form the coefficient of 
1/Pg (and 0{1/(5 G )) in eq. (|T3]). To understand better the behaviour in eq. (fH|), let 
us discuss observables simpler than ric- 

Consider a typical composite operator, such as (0*0). For (0*0), one can make a 
perturbative computation to find out what happens in the continuum limit. It turns 
out the there is a linear (1-loop) and logarithmic (2-loop) divergence. The finite MS- 
scheme continuum result (0*0(e|)) cont is 



(0*0(ej)) con t 1 . . 3.1759115^ 1 

^ - 2^(0 0)latt ^ 



log(6/3 G ) + 0.668 



(15) 



where "latt" refers to the normalisation of the field in the lattice action (^). The second 
term on the RHS is the linear, and the third the logarithmic divergence. The value 
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Figure 5: The distributions of the volume average of |^ GXj g G | for values of y just below 
and above the critical point y c at x = 0.0463. Precisely at y c , the two peaks would 
appear simultaneously. 

of (0*0)iatt measured on the lattice is thus not in itself a physical quantity: only its 
changes (such as the discontinuity across a first order transition) are, since in them the 
divergent parts cancel. 

We now expect that a similar thing happens for rtc- The difference is that nc is 
a more complicated and non-local quantity. The regularization sensitive part is not 
easily computable in perturbation theory, since even the integral 

N f S ce = J Vcf) \n c \ exp(- f d 3 x\d l( p\ 2 ) (16) 

is not Gaussian, due to the non-polynomial expression of \nc\- A numerical evaluation 
of N^ ee (a mass term has to be included on the lattice to kill a zero mode) gives a 
result close to the fitted values of f{n) in eq. flTBD (and favours c = 0). 

To demonstrate that the changes in N/3 gX /3 g are physical, the distributions of Np GX p G 
at x = 0.0463, (3g = 2, 4, 8 around the first order phase transition are shown in Fig. |5|. 
It is seen that the two-peak structure (in particular, the distance between the peaks) 
indeed remains the same within statistical errors when (3g is varied, when (3g is large 
enough (/5g = 4, 8). Thus the two-peak structure is physical and has a finite continuum 
limit, while the location of the structure on the \nj3 GX p G \-axis is unphysical and dom- 
inated by UV-effects: both peaks move to the right when (3g increases (the location 
of the y = 0.14 peak is shown in Fig. |j). Note that (3g = 2 is not yet in the scaling 
regime, and thus the distance between the peaks is different from that at (3g = 4, 8. 
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For x = 2, there is only a single peak which moves continuously to larger values as 
y is increased, see Fig. 0. 



5 Conclusions 

In conclusion, we have given a gauge-invariant definition for a vortex passing through 
a loop on a lattice, measured the corresponding total vortex density, and discussed 
its extrapolation to the continuum limit. We have pointed out that to approach a 
meaningful continuum limit, one must keep the size of the loop fixed in physical units. 
We have found that the total vortex density behaves as a pseudo order parameter, 
analogously to (4>*4>)~- the absolute value is always non-zero and is dominated by regu- 
larization effects near the continuum limit. Thus only the changes of the total vortex 
density with respect to the continuum parameters are physically meaningful. In the 
type I regime, the total vortex density displays a discontinuity, whereas in the type II 
regime, it behaves smoothly as the phase transition is crossed (Fig. |2|). 

The system possesses also (non-local) observables which behave analogously to true 
order parameters. One such is the photon mass, which vanishes exactly in the sym- 
metric phase. It has been suggested that another such quantity might be the density of 



long vortices passing through the whole lattice |23, [2(J . In contrast to the photon mass, 
this quantity should vanish in the broken phase and remain non-zero in the symmetric 
phase. The measurement of the density of long vortices is in progress fl9|j . 

Finally, it would be interesting to study the spatial distribution of vortices. This 
can be done by measuring correlators of "vortex density operators" separated by 
a distance r. One can define several correlators, depending on the relative orientations 
of the loops used in the nc's. These quantities would give realistic initial conditions 
for simulations of the time evolution of vortex networks. 
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